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ABSTRACT 



There is a long-standing controversy about the convergence of the dipole moment of the 
galaxy angular distribution (the so-called clustering dipole). Is the dipole convergent at all, 
and if so, what is the scale of the convergence? We study the growth of the clustering dipole 
of galaxies as a function of the limiting flux of the sample from the Two Micron All Sky 
Survey (2MASS). Contrary to some earlier claims, we flnd that the dipole does not converge 
before the completeness limit of the 2MASS Extended Source Catalog, i.e. up to 13.5 mag in 
the near-infrared Kg band (equivalent to an effective distance of 300 Mpc/h). We compare 
the observed growth of the dipole with the theoretically expected, conditional one (i.e., 
given the velocity of the Local Group relative to the CMB), for the ACDM power spectrum 
and cosmological parameters constrained by WMAP. The observed growth turns out to 
be within la confldence level of its theoretical counterpart once the proper observational 
window of the 2MASS flux-limited catalog is included. For a contrast, if the adopted 
window is a top-hat, then the predicted dipole grows significantly faster and converges 
(within the errors) to its final value for a distance of about 300 Mpc/h. By comparing the 
observational windows, we show that for a given flux limit and a corresponding distance 
limit, the 2MASS flux-weighted window passes less large-scale signal than the top-hat 
one. We conclude that the growth of the 2MASS dipole for effective distances greater 
than 200 Mpc/h is only apparent. On the other hand, for a distance of 80 Mpc/h (mean 
depth of the 2MASS Redshift Survey) and the ACDM power spectrum, the true dipole is 
expected to reach only ~ 80% of its final value. Eventually, since for the window function 
of 2MASS the predicted growth is consistent with the observed one, we can compare the 
two to evaluate /3 = Q'^^^/b. The result is ^ = 0.38 ± 0.04, which leads to an estimate of 
the density parameter Qm = 0.20 ± 0.08. 

Subject headings: cosmology: observations — cosmology: theory — (cosmology:) large-scale 
structure of universe — surveys — methods: data analysis — methods: statistical 
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Introduction 



The currently accepted paradigm of large-scale structure formation in the Universe is the 
gravitational instability. Within this framework, inhomogeneities in matter distribution induce 
gravitational accelerations, which result in galaxies having peculiar velocities that add to the 
Hubble flow. These velocities in turn enhance the growth of the inhomogeneities, causing strong 
coupling between cosmic velocity field and large-scale matter distribution. In perturbation theory of 
Friedman-Lemaitre models, in the linear regim e, peculiar ve locities and accelerations are aligned and 



proportional to each other at every point (e.g. 



Peebles 



1980) ■ More importantly, the proportionality 



coefficient of this relation is a simple function of the cosmological parameter of non-relativistic matter 
density, Q^, and practically does not depend on the cosmological constant (nor other forms of dark 
energy). Additionally, as the relation holds also for the motion of the whole Local Group of galaxies 
(LG) through the Universe, measurements of its peculiar velocity and acceleration may be used as a 
tool to dynamically constrain the parameter. 

The peculiar velocity of t he LG is known from the observed dipole anisotropy of the cosmic 
microwave background (CMB, 



to the barycenter of the LG fICourteau fc van den Bergh 



Hinshaw et al 



20091) ■ interprete d ClS db kinematic effect, and reduced 



19991 ). It equals to vcmb = 622 ± 35 km/s 



and points in the direction (/, b) = (272° ± 3°, 28° ± 5°) in Galactic coordinates. This kinematic 
interpretation is strongly supported by the observed alignment of the CMB dipole with the direction 
of the peculiar acceleration of the LG, although the latter is much more difficult to estimate. It 
requires knowledge of mass distribution in our cosmic neighborhood, and its determination had not 
been possible until deep all-sky galaxy catalogs became available. For that reaso n, the first attemp ts 



to measure the accele ration of the LG were made not earlier than 30 years ago flYahil et al 



1980 



Davis fc Huchra 



19821). 



Using an all-sky catalog, such measurement can be made under the assumption that visible 
(luminous) matter is a good tracer of the underlying density field. The general procedure is to 
estimate the so-called clustering dipole of a galaxy survey and infer the acceleration of the LG. 
However, such inference requires several conditions to be met. First, the survey should cover the 
whole sky; second, the observational proxy of the gravitational force (most often the flux of the galaxy 
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in the photometric band of the survey) should have known properties; and last but not least, the 
survey should be deep enough for the dipole to converge to the final value that we want to find. As 
usually one or more of these assumptions do not hold, the clustering dipole is a biased estimator of 
the acceleration and the estimation of the latter from the former may be done only if the mentioned 
effects are properly accounted for. 

In the past, many different datasets have been applied to calculate the clustering dipole. 
Generally speaking, there is no consistency on the amplitude, the scale of convergence of the dipole, 
and even on the convergence itself. The pioneering works used the revised Shapley-Ames catalog 



f Yahil et al. 



1980l ) and the CfA catalog (iDavis Sz Huchra]ll982l ). A great advancement came with the 



launch of the far-infrared IRAS satellite and catalogs obtained thanks to this mission. The LG dip ole 



from IRAS was studied first from two-dimensional d ata only (lYahil et al. 



1986 



1992 



Harmon et al. 



1987 



Schmoldt et al. 



1999 



Villumsen fc Strauss 



Rowan- Robinson et al.i 120001: 



20061 ) and with optical da ta added f|Lahav et al. 



data only were also used ( Lahav 



1995 



Kocevski fc Ebelins 



1987 



Hudso 



1 



1986 



Meiksin fc Davis 



19871). then with redshifts inc 



1988 



D' Mellow et al 



Kaiser fc Lahav 



2004 



uded ( 



Strauss et al. 



Basilakos fc Plionis 



19891 ). Sar nples with optic al 



19931 ). as well as galaxy clusters (IBrunozzi et al 



20061 ) . Among the most recent analyses one finds those directly related 



to the study presente d here, which used t he data from t 



Skrutskie et al. 



Catalog (XSC, 



20061: see also here under) . 



Mailer et al 



he Tw o Micron All Sky Survey (2MASS, 



Jarrett et al 



(120031 ) used the 2MASS Extended Source 



20001), concluded convergence of the clustering dipole from fiux data 



only and used it to calculate the average ma ss-to-light ratio in the near-infrared Kg band and 



to estimate the linear biasing parameter 6^. 



Erdogdu et al 



(120061 ) studied the acceleration of 



the LG from the 2 MASS Redshi 
measured redshifts, Huchra et al. 



' t Sur vey (2MRS, an all-sky subset of the 2MASS catalog with 



20051 ) and estimated a combination of the mass dens i ty flm and 



Lavauxetal 



2010|) used 



luminosity bias parameters, the so-called (3 parameter. In a recent work, 
an orbit-reconstruction algorithm to generate the peculiar velocity field for the 2MRS, extended it to 
larger radii, and observed no convergence up to at least 120 Mpc/h. 

Since both the gravitational force and the fiux of a galaxy are inversely proportional to the 
distance squared, the clustering dipole can be calculated by assuming some average value of the 
mass-to-light ratio and summing only fiux 'vectors' of galaxies (with positions on the sky used as 
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angular coordinates). Such a dipole does not use distance nor redshift information of particular 
objects, which hinders the estimation of the depth up to which the dipole is measured and can bias 
the estimation of LG acceleration from it. This drawback is partially removed if a galaxy redshift 
survey (with redshifts of galaxies used as proxies of their distances) is used for the analysis; however, 
up to date the densest and deepest all-sky redshift surveys do not reach farther than to ~ 100 Mpc/h 
and contain no more than s everal dozen thousand sources. Owing to som e recent claims of large-scale 



flows in the local Universe (IWatkins et al 



2009 



Kashlinsky et al. 



20081 ). it is important to estimate 



the dipole from as big catalogs as possible, even if they do not contain redshift measurements. We 
thus have decided to trade the advantages of redshift measurements for a huge number of galaxies and 
a much greater depth of a photometric-only catalog. The angular dipole thus o btained is ad ditionally 



free from any so-called redshift distortions, in particular from the rocket effect ([Kaiser 



19871 ). 



In this Paper, we focus mainly on the issue of the convergence of the dipole, by analyzing the 
growth of the dipole with increased depth of the sample. For that purpose w e use the data from 



Mailer et al. 



toom . One of the goals of 



the 2MASS XSC, similarly to what has been presented in 
this Paper is to expand and refine the latter work. We do it by modifying the criteria of removing 
Galactic objects from the sample and masking and filling the Galactic plane and bulge (the Zone 
of Avoidance, ZoA), but keeping in the ZoA the brightest and closest galaxies of big influence for 
the local motion. We then analyze the growth of the dipole and check how its direction changes on 
the sky. Next, we compare the observed growth of the dipole with theoretical expectations, using 
the condition of kn own peculiar velocity of the Local Group, in a similar manner as was done in 



Lavaux et al. 



(I2OIOI ). From this comparison we evaluate the combination of the cosmological matter 
density parameter and the linear bias (the '/?' parameter) and, by using externally provided value of 
bias, constrain the parameter f^m. 

The paper is organized as follows. In Section [2l we present the theoretical framework for the 
calculation of the clustering dipole. Section [3] contains details on how we handled the 2MASS data. 
Next, Section m focuses on the growth of the clustering dipole. Section [5] focuses on the position of 
the 2MASS dipole on the sky. In subsequent Section Owe discuss the results: first, in Subsection 16.11 
we present the theoretical basis to analyze the growth for known LG velocity; Subsection 16.21 focuses 
on the observational window of our measurement; finally in Subsection 16.31 we compare our results 
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with theoretical expectations and use this comparison to estimate the /3 parameter from the growth 
of the dipole alone. We shortly summarize and conclude in Section [TJ where we also provide the 
obtained value of the f2m parameter. 



2. The clustering dipole 



The gravitational instability scenario of large-scale structure formation relates peculiar velocities 
of galaxies with their p eculiar g r avitat ional accelerations. In linear theory, this relation has a 



particularly simple form 



PeeblesI (11980f ): 



V 



9 



2f{Q, 



9 



Here, Hq = 100 h km/s/Mpc is the Hubble constant, is the current value of the cosmological 
density parameter of non-relativistic matter, /(^m) = (dhiD/dlnz) \z=q (with D being the growth 
factor) and pb is the mean matter density of the background. Within ACDM mod els with a 



cosmological constant, the growth p arameter f is very well fitted by f{^m) — t jLinder 



and is virtually independent of A, cf. 
given by the integral: 



Lahav et al. 



20051 ) 



(jl99lh . The acceleration vector at a position 'r' is 

(2) 



where Sm{r) = [pm(''") — Pb] / Pb is the density contrast of non-relativistic matter at the point r. 
However, as what we observe are galaxies, we have to assume some relation between their density 
field and that of matter. This is usually done via the linear biasing paradigm: 5^ = b6^. This biasing 
scheme, valid in linear theory that we use in our whole analysis, neglects the stochasticity, as well as 
possible scale- and galaxy-type dependence in the relation between the two density fields . For more 
details on possible non-linear biasing, see for example the review by iLahav fc Sutd (120041 ). 



Including the biasing relation into Eq. ([2]) and using the fact that for a spherical survey 
/ \r'-r\'-^ d^?"' = 0, we get the following expression for the peculiar acceleration: 



9ir) = ^ 



dV . 



(3) 



The biasing parameter b is usually combined with the factor /(i^m) into the parameter /3 = /(f2m)/6. 
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Comparing Eqs. ([T]) and ([3]), we get the proportionality valid in linear theory: 

v(xPg. (4) 

This relation could in principle be applied directly to the motion of the Local Group of galaxies 
through the Universe. Consequently, comparison of the peculiar velocity and acceleration of the LG 
serves as a tool to estimate the (3 parameter. Independent knowledge of biasing allows to estimate 
cosmological density VL^. 

In reality, however, we do not observe continuous galaxy density field, but rather discrete 



objects, even if in a 



yery 



arge number. In the following derivation, which can be found e.g. in 
Villumsen &: StraussI (119871 ). we model galaxies as point sources: Pg{r) = ^ . Mj Sd^t — rj), where 6^ 
is Dirac's delta; Mj and are respectively the mass and the position of the i-th galaxy. Putting the 
coordinate system at r = (barycenter of the LG) we obtain the acceleration of the LG as a sum of 
the force contributions from all sources in the Universe: 



This Newtonian formulsLj is still not useful for calculations based on observational data, as masses of 
individual galaxies are usually known very poorly, if ever. However, if the i-th galaxy has an intrinsic 
luminosity Lj, we can write 

i ' i 

where Si = Li/Anrf is the flux received from the z-th object. This relation means that if we know the 
behavior of the mass-to-light ratio in the band(s) of the survey, we can estimate the acceleration of the 
LG from a two-dimensional catalog, i.e. one containing astro- and photometric data only (positions 
and fluxes). Furthermore, if the mean mass-to-light ratio is a universal constant, T = (M/L), we 



^The Newtonian limit can be applied as our whole analysis concerns distances well below the Hubble 
radius rn = c/Ho = 3Gpc//i. 



finally ge 

i 

In some applications, including the present one, it is more convenient to work in terms of matter and 
luminosity densities. This is especially the case when the luminosity density, j, is known for a given 
band, rather than the mass-to-light ratio. We have 



L/ J SttGj 
which gives 

« = ^ES.^.. (9) 

The luminosity density j for a particular band of the survey can be calculated for example from the 



luminosity function ^{L) of galaxies in this band (e.g. |Peebleslll993l ) 



J 



oo 

J L$(L)dL . (10) 



Note that using the Relation ([9]) in Eq. ([T]), we get the linear-theory velocity measured from the flux 
dipole as 

^ = /3 — V5,f, = /3^, (11) 

where g denotes the scaled acceleration of the LG, in units of velocity. The term SiVi is the flux 
dipole moment of all sources down to the zero flux over the whole sky. The universal luminosity 
density j, measured from a fair sample of galaxies in the given band with known apparent luminosities 
and redshifts, is proportional to Ho, which means that the overall result does not depend on the 
Hubble constant. 

Realistic galaxy catalogs will never reach down to zero flux, irrespectively of the used wavelength. 
On the contrary, surveys are usually flux-limited, which means that the number of observed sources, 
A^, is finite. For that reason, in the following we will denote the flux dipole of a finite sample as d: 

3 Qm H? 



2bj 



(12) 



^Note that as we are interested in the motion of the Local Group as an entire system, the galaxies 
of the LG should not be included in the summation. 
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Note that the clustering dipole calculated for a finite, fiux-limited sample may be a biased estimator 
of the peculiar acceleration of the Local Group, Eg. This can be overcome by extrapolating the 



measured dipole to zero fiux (IVillumsen fc Strauss 



19871 ). The situation gets easier in the case of 



galaxy redshift surveys. We can then use the redshifts as distance estima tors and weight galaxies with 



the inverse of the selection function of the survey (e.g. lYahil et al.lll99ll ). in order to mimic an ideal, 
volume-limited catalog. However, despite an outstanding advancement in surveying the cosmos in 
recent years, the deepest and densest all-sky redshift survey, 2MRS, contains only ~ 43,000 galaxies 
and has a median depth of merely ~ 100 Mpc/h {z^ ~ 0.04, Huchra et al. 2011, in preparatioio) . On 
the other hand, the 'parent' catalog of this survey, namely 2MASS XSC, reaches 3 times deeper and 
includes o ver 20 times more galaxies. Its redshift coverage , when matched with other surveys, such 



as SDSS (Aihara et al. 



201 ih or 6dFGS flJones et al 



depth (for the most recent compilation see 



20091). is non- uniform both on the sky and in 



Lavaux fc Hudson 



20111). For the purpose of our work, we 



have decided to sacrifice the advantages of weighting galaxies, possible for redshift surveys, obtaining 
instead an overwhelmingly greater number of sources and unprecedented depth of the survey with 
photometric data only. An additional motivation of using the dipole f[T^ . constructed only with the 
us of fluxes and angular positions of individual gal axies, is the fact that it is free of any redshift-space 



distortions, and in particular of the rocket effect (IKaiser 



19871 ). The latter consists in the fact that 



the peculiar acceleration of the LG calculated using redshifts instead of real distances will differ from 
the actual LG acceleration due to a spurious contribution from the galaxies that are in the direction 
of the LG motion. Here, we do not use distances measured in redshift- nor in real space, and the only 
possible effect of that kind would be the anisotropy modulation in the distant gala xy distribution 



due to the aberration effect, which is however completely negligible for our sample (jitoh et al. 



2010h . 



The only stage at which the Kaiser effect comes into play is in the measurement of the luminosity 



function $(L) and consequently 
where $(-L) is measured, see e.g. 



:he luminosity den sity j. This is addressed in the relevant papers 



Jones et al 



(120061). 



The main interest of this paper is to examine the dipole f|T2l) as a function of increasing depth of 
the 2MASS galaxy sample, as described in detail in the following chapters. This analysis will allow us 



^For current status see http://tdc-www.cfa.harvard.edu/2mrs/. 
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to constrain the /3 parameter from this dipole and consequently to measure the fim parameter. Our 
method however is not to directly compare the peculiar velocity and acceleration of the Local Group; 
instead, we use the observed growth of the dipole to obtain these constraints. 



3. Data preparation 



The Two Micron All Sky Survey (2MASS, ISkrutskie et al.l 120061 ) is the first near-infrared 



survey of the whole sky (covering 99.998% of the celestial sphere), and was performed in the period 
1997-2001 in the J(1.25/im), if (1.65 /xm) and Kg {2.16 fim) bands, with the use of twin 1.3-m 
ground-based telescopes. All the data from the survey are available through the NASA/IPAC 
Infrared Science Archived The main outcome of this project are two photometric catalogs: of point 
sources (PSC), containing about 471 million objects, and of extended ones ( XSC), with in ore than 1.6 



million objects, mainly galaxies (> 98%) and some diffuse Galactic sources (jJarrettll2004l ). The XSC, 
which was used for the purpose of our analysis, is complete for sources brighter than Kg ~ 13.5 mag 
(~ 2.7 mJy) and resolved diameters larger than ~ 10 - 15". The near-infrared flux is particularly 
useful for the purpose of large-scale structure studies as it samples the old stellar populati on, anc 



hence the bulk of stellar mass, and it is minimally affected by dust in the Galactic plane (jJarrett 



2004J ). An additional advantage of using 2MASS data, especially in the context of calculating the flux 
dipole, for which apparent magnitudes are used, is the global photometric uniformity of the catalog, 
which was enforced by nightly photometric calibration to an extensive set of standard star fields. 
On the other hand, as any survey, 2MASS is not perfect. It is biased against optically blue and low 
surface brightness galaxies, such as dwarfs, but sensitive to the early type, bulge-dominated ones. As 
the former have very small luminosities and masses, their possible underrepresentation in the catalog 
should not influence significantly our results. 

The 2MASS photometry offers several types of 'magnitudes' for extended objects, depending 
on the type of aperture used etc. Throughout the whole analysis we use the 20mag/sq." isophotal 
fiducial elliptical aperture magnitudes, which are defined as magnitudes inside the elliptical isophote 



^http : / /irsa . ipac . caltech . edu/Missions/2mass . html 
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corresponding to a surface brightness of /^band = 20 mag/sq.". We prefer those to the Kron ones as the 
latter use large and noisy apertures, prone to contamination, resulting in sys tematic overestimatioii . 
Our choice is additionally supported by the considerations in the appendix of iKochanek et al.l (120011 ) . 
However, we must remember to correct the values used by adding an offset of A = —0.2 mag when 
converting to flux (cf. S ection HI) , in order to compensate for the flux lost outside the aperture 



(typically ~ 10% - 20%, 



Jarrett et al. 



20031 ). We have checked that this offset is roughly equal to the 



one between isophotal fiducial elliptical aperture magnitudes and the 'total' ones, obtained from fit 
extrapolation. The magnitude correction by a constant factor certainly introduces some scatter in 
total flux estimates, as in reality it may depend on galaxy morphology. The latter is very hard to 
constrain from 2MASS data, we will thus treat this scatter as a systematic effect that needs to be 
included in the error budget. Our tests show that this error is of the order of a couple percent. 



In order to prepare the data for our purposes, we hay e proceeded as follows. First of all, we 



applied the extinction correction from 



Schlegel et al. 



( 1998 ). by calling the procedure dust_getva: 



I 



for galactic coordinates of each of the objects. The procedure yielded values of E{B — V), which 
were subtracted from the original magnitudes with appropriate multiplicative factors Ry taken from 



Cardelli et al. 



(119891 ): 0.902 for J, 0.576 for H and 0.367 for K. W e performed the subtraction for 



objects with \b\ > 5° due to the statement of 



Schlegel et al 



mm that for |6| < 5° the predicted 



reddenings should not be trusted (the sources in the ZoA were not included in our catalog apart from 
several brightest; see end of this Section). Moreover, for some minor parts of the sky, the extinction 
correction gave unreliably high reddenings, which resulted in some objects becoming unrealistically 
bright (with negative magnitudes) and eventually deleted. At this stage, we have also identified and 
removed the following objects (with some found in more than one category): 

- artifacts: flag cc_f lg=a in the 2MASS XSC (122 objects); 

- sources with NULL or unreliable K magnitudes (as described above) (718 objects); 



- non-extended sources: flag vc=2 in the 2MASS XSC (7383 objects); 



^http : / /www . astro . princeton . edu/~schlegel/ dust/ dustpub/ CodeC/README . C 



- Local Group galaxies, taken from the list of 



Lee fc Leel fl2008f ) (31 objects) 



,0 



Milky Way sources, identified by cross-correlation with a list of such 4454 objects, separately 
identified earher in the 2MASS XSC. 



The 2MASS Extended Source Catalog contains mainly galaxies; however, it is also comprised of 
Milky Way entities, such as stellar clusters, planetary nebulae, HIT regions, young stellar objects and 
so on. In order to keep our analysis reliable, these objects had to be removed from the catalog. This 
was partially done for the 4454 sources mentioned above. However, owing to the size of the catalog, 
any further 'manual' procedure of Galactic object removal was impossible and only a method based 
on some general properties could be applie d. A useful one in t his regard is the color, i.e. difference of 



magnitudes in two bands. In their analysis. 



Mailer et al. 



(120031 ) made a cross-correlation with galaxies 



spectroscopically confirmed from the Sloan Digital Sky Survey (SDSS) and excluded extended sources 
brighter than = 12 mag with colors J — K < 0.75 or J — K > 1.4; at fainter magnitudes only those 
objects with J — K < 0.5 were removed. We have checked these conditions by examining the J — K 
distribution of galaxies in the 2MASS Redshift Survey (2MRS), which is a complete subset of the 



2MASS XSC, containing all the galaxies with K < 11.75 mag, with measure d redshifts 



43,000 objects). Using the data from the Extragalactic Distance Databasqj ( iTully et al. 



more than 



20091), we 



have found that in deed galaxies are cl ustered around J — K = 1; however, we have decided to alter 



the limits given by 



Mailer et al. 



3). 



(120031 ). Analyzing additionally the distribution of XSC objects with 



K < 13.5 mag and |6| > 15°, among which there are mainly non-Milky Way sources, apart from some 
molecular clouds (Kaluzny, private communic ation), we ha ve decided to keep in our catalog those 



objects that have 0.6 < J — K < 2.0 (cf. also 



Jarrett 



2000) • We us e this criterion for all sources 



independently of magnitude, as we think that a differentiation as in 



Mailer et al 



(l2003h could lead to 



a bias in the sample. An additional eyeball verification of 100 brightest objects which pass this filter 



^Not all the objects from the Local Group were found in the database. These were some dwarf 
galaxies of low mass and near-IR luminosity, hidden behind the Galaxy or with surface brightness 
below the threshold of 2MASS. 

^http : //edd. if a.hawaii . edu; the complete 2MRS catalog is currently downloadable from 
http : / /tdc-www . cf a . harvard . edu/ 2mrs/. 
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off the Galactic plane confirms that indeed all of them are galaxies and that only one galaxy with an 
extreme value of J — is removed by this procedure up to i^^ ~ 7.5 mag. 

An important issue in the calculation of the clustering dipole is the Zone of Avoidance (ZoA), 
i.e. the region of the sky with small Galactic latitudes b, which obscures galaxies behind the Galactic 
plane and bulge. Although the Galactic extinction is much low er in the near infrar ed than in visible 



bands fiCardelli et al.iil989i ) and this applies equally to the ZoA (iJarrett et al. 



2OO0I), the 2MASS XSC 



is sti ll incomplete near the Galactic e quator, mainly due to high stellar density in this region of the 



sky (IKraan-Korteweg fc Jarrett 



maps of 



Schlegel et al. 



1998 



2OO5I ) . For that reason, and owing to inapplicability of the extinction 



ane and bulge in the 



Erdogdu et al. 



mm, 



for |6| < 5°, we have masked out the Galactic p 
following way. For the shape of the mask we have chosen the one proposed in 
i.e. we have skipped all the objects with \b\ < 5° (plane) and \b\ < 10° for I < 30° or / > 330° (bulge). 
Then we have filled the resultant gap by cloning the adjacent strips, with mirror-like reflections: 
for instance, objects with 10° < b < 20° were copied to the bulge by assigning fenew := 20° — 6oid 
and keeping other parameters unchanged (such as the longitude / and magnitudes). An analogous 



advantage over random filling (considered both in 



for the Ga 


actic 


Mailer et al. 


2003 



actic plan e. Such cloning has the 



and 



Erdogdu et al. 



20061 ) that it 



extends the structures from above and below the ZoA; moreover, the only artificial discontinuity of 
the galaxy distribution created in this procedure is at the Galactic equator and at the edges of the 
box masking the bulge. We have also tried other masks and methods of filling the ZoA, and found 
no special importance for the results of the analysis presented here. This will be shortly addressed in 
Section HI 

Once the ZoA has been masked and filled, we have added to our catalog several galaxies that 
were not present in the 2MASS XSC but could be found in the 2MASS Large Galaxy Atlas (LGA, 



Jarrett et al 



20031 ). This atlaqj contains the ~ 600 largest galaxies as seen in the near-infrared, of 
which some 50 are not present in XSC or are located in the ZoA, \b\ < 5° (17 sources). Among the 
latter, three are of particular importance for the Local Group motion, namely Maffei 1 & 2 and 
Circinus. We will discuss their influence on our results later in the text. Note that this addition does 



^Accessible through IRSA at http://irsa.ipac.caltech.edu/applications/2MASS/LGA/ 
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not spoil the photometric uniformity of the resulting sample because all the galaxies from the LGA 
present also in the XSC were assigned the magnitudes from the former catalog when the latter one 
was constructed. 



In case of those LGA galaxies that were present in the ZoA, the 



Schlegel et al. 



are known to overestimate the Galactic extinction by roughly 15% (e.g. 



Schroder et al. 



1998h m aps 



20071 ). we 



have thus decreased the E{B — V) by that amou nt there. The ex c eption s are Maffeis 1 & 2, for 



which we used exact values of extinction, given in 



Fingerhut et al 



(120071 ). as well as Circinus with 



E{B — V) = 0.677 (For et al. 2011, in preparation). Note however that apart from those three 
galaxies, which have extinction-corrected magnitudes below K = 5 mag, all the remaining ones added 
from LGA are much fainter, at least by 2 mag, and possible misestimation of their extinction does 
not largely influence our analysis. 



4. Growth of the 2MASS dipole 

In this section we will use the data prepared as described in Section [3] to calculate the clustering 
dipole of the galaxies from the 2MASS XSC and analyze its growth. For that purpose we change 
units of d into km/s, defining the scaled dipole as 

d=:rp^d = ^f2s.n, (13) 

SHqII^ bj ^ 
where we have used the Formula (fT2l) for the flux dipole of the survey. 

The flux of each galaxy is calculated from its magnitude as 

Si = So lO-o-*'"' , (14) 

where 5*0 is the flux for a 0-magnitude object. As was already stated, we consider magnitudes in the 
Kg band, which was the main ('target') band of the 2MASS survey (for simplicity of notation, we 
sometimes skip the 's' subscript). Then Eq. (IT^ takes on the form: 

d=J^ S{0 mag) y io-o.4(i^.+zPO)r^ ^ (15) 
bjK Y 

where gfOmag) = 1.122 x lO^^^Wm-^ (±1.685%) and the zero point offset ZPO = 0.017 ± 0.005 



(ICohen et al. 



20031 ). The Ki magnitudes in Eq. (ITSl) include also a negative offset of A = —0.2 mag 
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added following iKochanek et al.l (l200ll ) due to the underestimation of total fluxes by the isophotal 
magnitudes in the 2MASS XSC The quantity jx is the luminosity density in the K band. It is 
obtained for example from the integral (fTOj) using the luminosity function in this band. The value 



of ?V h as been est i mate d 



(I2QO1I); 



Bell et al. 



(120031 ) 



by many author s in the recent decade: iKochanek et al 



Eke et al. 



(120051 ) 



Jones et al 



one we consider the most reliable, calculated by 



1120061) 



Jones et al. 



Smith et al 



(1200 ll ) 



Cole et al. 



(I2OO9I ). We use here the 



(120061 ) from the 



of more than 60,000 galaxies in the 6-degree Field Galaxy S urvey (6dFGS, 



jx = (5.9 ± 0.6) X lO^/iLf Mpc-^ with Lf ~ 6.8 x lO^^W mieke et al.lboosh pl Note that as the 



luminosity funct ion 



Jones et al 



20091): 



luminosity density depends on the Hubble constant, the h factor cancels out and Eq. ([T^ may be 
rewritten as 



C 



OAK,' 



(16) 



where C ~ 2620 km/s for parameter values as given above. The limits of the summation are now the 
minimum and maximum K magnitudes of galaxies in the sample. The lower limit in Eq. (fT6|) is the 
magnitude of the brightest object in the catalog (excluding LG galaxies) and the upper one will be 
increased, as will be discussed later in the text (Section Hj). To retain the reliability of the analysis. 



we must note that the results remain trustworthy only up to th e comp 



the sample becomes incomplete for objects with K > 13.5 mag (jjarrett 



eteness limit of the catalog: 



2004 ). 



The growth of the dipole was calculated by incrementing -ft^max in the sum given in Eq. ([T^ . 
Results are illustrated in Figures [1] and [2l Figure [1] shows the growth of the dipole as a function 
of the number of galaxies used for the calculation, together with Galactic Cartesian components 
of the acceleration [g^ points towards the Galactic center, I = 0° , b = 0°; Qy is perpendicular to 
it in the Galactic plane ; is perpendicular to the Galactic plane). The similarity of this plot to 



figure 1 of 



Mailer et al. 



(I2OO3I ) is intentional a nd allows us t o dra w the conclusion that although our 



data analysis was slightly different to that of 



Mailer et al 



(I2OO3I ) (altered shape and filling of the 



ZoA, different removal of Galactic objects), the two approaches give qualitatively the same results. 
However, we want to stress that such a presentation of the data may be misleading. Owing to local 



^Lq is the Solar luminosity in the K band, calculated based on table 7 of [ Rieke et al. 
http : //www . iop . org/EJ/article/1538-3881/135/6/2245/aj 271287_mrt7 . txt. 



(120081 ) 
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8 



number of galaxies (x10 ) 



Fig. 1. — Growth of the 2MASS clustering dipole as a function of the number of galaxies used for 
the calculation (bottom axis), ordered by their Kg magnitudes (top axis). The thick red line is 
the amplitude of the dipole; the thin dotted and dashed lines (green, blue and magenta) are the 
Cartesian components (in Galactic coordinates). Two vertical lines illustrate the limits of the 2MRS 
Ks < 11.25 mag sample and completeness of the 2MASS XSC {Kg = 13.5 mag). Such data presentation 
could suggest that the 2MASS dipole has converged within sample limits, which is not the case. 
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Fig. 2. — Growth of the 2MASS clustering dipole as a function of increased maximum Kg magnitude 
of the subsample (bottom axis). Corresponding effective distance is given at the top. The thick red 
hne is the amphtude of the dipole; the thin dotted and dashed lines (green, blue and magenta) are 
the Cartesian components (in Galactic coordinates). Two vertical lines illustrate the limits of the 
2MRS Ks < 11.25 mag sample and completeness of the 2MASS XSC [Kg = 13.5 mag). The lack of 
convergence of the dipole is clearly visible. Note also the steady growth of the amplitude as well as of 
the Galactic y component of the dipole. 
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inhomogeneities, a linear scale in the number of galaxies on the x-axis does not have to be related 
in a simple manner to a scale expressed in cut-ofF magnitudes nor in distances and 'compresses' 
the left-hand side of the plot, while 'stretching' the right-h and one. This results in a suggestion of 



convergence of the dipole, which was actually concluded by iMaller et al.l (120031 ). However, what we 
actually need in order to draw proper conclusions about the growth and possible convergence of the 
dipole, is a linear scale in distance on the abscissa. 

Neither distances, nor even redshifts (including photometric ones) are currently meas ured for 



the whol e 2MASS XSC (alt 



2004 and 



lough some attempts are being made regarding photo-z's, see 



Francis Sz Peacock 



Jarrett 



2010l ). We thus need to deduce effective distances of galaxies from their 
fluxes, with the use of the luminosity function (LF) in the K band. If all the galaxies had the same 
luminosity, say L^, the relation between the observed flux S and distance r would be straightforward: 
r = a/ L^./AttS. However, galaxies have different morphologies, masses and luminosities, and obviously 
their LF is not a Dirac's delta (on the contrary, it is very broad). Therefore an estimated distance 
of a galaxy with a given flux must have a scatter. Constructing the estimator, a first choice could 
be the conditional mean, i.e. the expectation value for r given 5*. Instead, we think that it is better 
to choose the conditional median for r^.^ (a median value of distance given the flux). We consider it 
being more adequate to our problem: the same number of galaxies with a given flux have distances 



smaller and greater than the median. 



Detailed calculations , presented i n Ap 



pendix 



Jones et al. 



K band LF as given by 

= -23.83 + 51og/i ± 0.03 and a 
equals to 



(see also 



Peebles 



(120061 ) ■ i.e. a Schechter function (jSchechterl Il976l ) with 



19931). show; that for the 



-1.16 ± 0.04, this effective distance for the magnitude K 



Teff ~ 0.59 X 10°-^^ Mpc/h . 



(17) 



This proxy of distance is used in Figure [21 which differs from Figure [H by a different scaling of the 
X-axis. The growth of the clustering dipole up to the completeness limit of the 2MASS XSC is now 
evident. Additionally, note that the growth has an essentially constant slope for > 150 Mpc/h, 
i.e. > 12 mae. 



An interesting feature is the behavior of the Galactic Cartesian components of the dipole. The x 
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and z ones are virtually constant for refr > 150 Mpc/Zi; however, the y component still grows even at 
the limit of the catalog, similarly as does the total amplitude. This could point to some systematic 
effect, related to masking and filling of the Zone of Avoidance. We have however checked that the 
same qualitative behavior of the three components is observed for different shapes of the mask and 
the way it is filled; what is more, the effect exists even if we calculate the dipole having removed from 
the catalog all the galaxies with |6| < 10° (leaving the resulting strip completely devoid of galaxies). 
We have also observed that adding to the sample the LGA galaxies that were not present in our 
catalog has virtually no influence on the amplitude of the dipole and only slightly changes the values 
of particular components (see also below). 



5. Misalignment angle 



The linear theory Relation ([T]) between peculiar velocity and acceleration of the Local Group 
predicts that the two vectors should be parallel. In reality, a non-zero misalignment angle between 



Crook et al. 


2007 


) or stochas 


ticity 


of the LG [\ 


3ernardeau et al. 


1999 



Cieciel^g et al 



2001 



Chodorowski 



of this angle are usually of the order of 10 - 20 degrees (IStrauss et al 



Mailer et al 



2003 



Erdogdu et al. 



et al. 



ocity and acceleration 



1992 



20081). Observed value s 



Schmoldt et al. 



1999 



20061 ) and we have obtained similar results, which confirms the 
validity of the linear approximation for the scales of interest. Figure [3] shows how the dipole direction 
on the sky changes as the sample depth increases, from 'min', i.e. 5 Mpc//i, to 'max', equal to 
310 Mpc//i. The black disc shows the CMB dipole. For comparison, we present also t wo ot her results 



Mailer et al. 



fl2003h fgr een 



Erdogdu et al. 



fl2006h 



of 2MASS data analysis: 2MASS dipole with 'cloning' the ZoA, from 
triangle labeled M03) and the flux weighted 2MRS dipole in the CMB frame of 
(magenta square, E06). We increase the Galactic longitude I from left to right for easier comparison 
with relevant figures in these two papers (respectively fig. 2 of M03 and fig. 7 of E06). 

In Figure H] we plot the misalignment angle with respect to the CMB dipole direction, as a 
function of growing depth of the sample. The two lines illustrate the effect of adding the LGA 
galaxies to the catalog, and in particular of the three bright galaxies located in the ZoA (Maffeis 1 & 
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Fig. 3. — Position of the 2MASS clustering dipole on the sky as a function of increased depth of the 
sample. The grid shows Galactic coordinates: longitude / (labels on bottom) and latitude b (left). 
The min/max labels refer to minimum and maximum depth of the subsample for which the dipole was 
calculated, respectively 5 Mpc/h and 310 Mpc/h. The black disc marked by 'CMB' is the direction 
of the peculiar velocity of the Local Group (/ = 272°, b = 28°). The green triangle labeled 'M03' 
shows the 2MASS clustering dipole of Mailer et al. (2003): / = 263°, b = 40° (with cloning in the 
ZoA). The magenta square marked by 'E06' is the direction of the flux-weighted 2MRS dipole in the 
CMB frame of Erdogdu et al. (2005): / = 245°, b = 39°. Blue labels indicate important structures in 
the local Universe: the Virgo Cluster {D ~ 17 Mpc), the Hydra Cluster {D ~ 47 Mpc) and the Leo 
Supercluster {z > 0.031). 
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Fig. 4. — Misalignment angle between the 2MASS clustering dipole and the peculiar velocity of the 
Local Group as a function of increased maximum Kg magnitude of the subsample (bottom axis). The 
corresponding effective distance is given at the top. Two vertical lines illustrate the limits of the 2MRS 
Ks < 11.25 mag sample and completeness of the 2MASS XSC {Kg = 13.5 mag). The solid and dashed 
lines are shown to illustrate the effect of adding the Large Galaxy Atlas to the sample. 
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2 and Circinus). A specific 'tug-of-war' between these structures, located on almost opposite sides of 
the Galactic plane, results in a general raise of the angle by ~ 1°. As can be also seen, the brightest 
(and presumably the closest) galaxies have the most influence on the value of the angle. Note also 
that already for a distance as small as 50 Mpc//i the misalignment angle reaches a value of ~ 20° 
(~ 0.35 rad), with a minimum of 16.8° for rgg = 305 Mpc/h (without the LGA galaxies). 

One can think of several reasons for this misalignment. The first possibility are influences from 
extremely larg e scale s. We cannot addres s this issue directly here, note only that the claims of 



Watkins et al 



fl2009h or 



Kashlinsky et al 



fl2008h of large-sca 



than e xpe cted in ACDM have 



(120 llh or 



Osborne et al. 



3een recently put in doubt by 



e bulk flows with va l ues much larger 



Nusser et al. 



fl2011bh 



Nusser fc Davis 



20101); the debate is thus far from being settled. The second option 



could be unsuspected s tructu res in the Zone of Avoidance. This issue was partly addressed in 



Bilicki &: Chodorowskil (120101 ). where the influence of one of such partially obscured objects - the 
Local Void (LV) - was analyzed. As was shown, the improper fllling of the ZoA in the LV region 
cannot bias the angle by more than 1°. Additionally, already in 2005 it was quite certain that 
thanks to surveys in various wavelengths (such as HI or X rays) there should b e no unknown nearby 



Fairall fc Lahav 



20051 ). Finally, such a big 



large-scale structures hidden behind the Milky Way (e.g. 
value of the misalignment angle can be due to improper accounting of very local influences; we have 
tested this possibility by adding the LGA galaxies to our sample. 

Note however that our statistical assumption that the flux is a good proxy of the gravitational 
force, which may be appropriate for the whole sample, can fail for individual galaxies. In particular, 
we have used a constant mass-to-light ratio for all galaxies, which is probably valid when averaged 
over many of them, especially since we use the Kg-hand luminosity, which is known to be a better 
tracer of stellar mass than optical and other NIR wavebands. Still, the M/L^ ratio should vary 
somewhat with morphology and luminosity. For galaxies with smallest magnitudes these variations 
may sign iflcantly affect thei r contribution to the direction of the clustering dipole. The same was 



noted by 



Mailer et al 



(120031 ). who also concluded that the 16° difference between their 2MASS dipole 
and the CMB one was caused by the brightest galaxies. We obtain a similar result here and interpret 
it as due to a speciflc non-linear effect, the shot noise, which translates into the fact that just a few 
sources with the largest fluxes have a large influence on the misalignment. 
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The high significance of the most luminous (in ter ms of the observed flux) galaxies for the 



direction of the dipole conflrms theoretical predictions of 



Chodorowski et al. 



(120081 ): the misalignment 



angle could be lowered if the observational window used was the optimal one, which would be the 



case if we removed those brigh test galaxies from the sample (as was also tried both by 



Mailer et al. 



2003 and 



Erdogdu et al. 



20061 ). We plan to address this issue in a forthcoming paper concerning the 



maximum-likelihood estimation of f3 from the 2MASS dipole. 



6. Discussion 



Our results showing that the clustering dipole of 2MASS galaxies still grows even at the 
completeness limit of the sample, taken at face value, are consistent with the results of some othe r 



autho rs, who used v arious datasets a r id me thods: redshift sur vey of IRAS galaxie s 



fll992l ): IRAS PSCz, 



selected clusters. 



Schmoldt et al. 



9991) and reanalysis by 



Kocevski fc Ebelind l\20m : reconstructed velocity fleld of 2MRS, 



Basilakos fc PhonisI ( 20061): X-rav 



Strauss et al 



Lavaux et al. 



(120101 ). On the other hand, t hey cont r adict claims of convergence at scales even as s mall as 60 



OOMpc//i: optical s ample of 



mm\ IRAS PSCz, 



HudsonI ( 



3 



'. lowan-Robinson et al. 



2MRS. lErdogdu et al 



Mailer et al 



1995 ) ; reds hift sample of Abell/AC O clusters. 



(I2OO0I ): IRAS FSCz and BTP, 



Brunozzi et al. 



D' Mellow et al. 



(12004f ) 



(I2OO6I ). Our analysis also suggests a different interpretation of the results of 



(I2OO3I ) — data presentation as in our Figure [2] instead of flgure 1 therein would possibly 



point to similar lack of convergence . On t he other hand, in order to ba able to directly compare our 



results with those of 



Erdogdu et al. 



( 120061 ). we would have to apply the same weighting as was done 
for the 2MRS sample, namely by the inverse of the flux-weighted selection function. We are unable 
to do it, not knowing distances nor redshifts for the whole sample. 

Apart from galaxy weighting, the discrepancies between the above listed results most probably 
stem from the different nature of catalogs and methods used for the calculation, and in particular 
may be due to distinct observational windows. Such a window for a given survey describes the sample: 
it may be interpreted as a fllter (in real or Fourier space) through which we observe the Universe. 
Knowledge of the observational windows, necessary to correctly confront results as those given above, 
is also essential if we want to make comparisons with theoretical expectations. 
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6.1. Theoretical framework 



We would now like to check if the behavior of the 2MASS flux dipole is consistent with the 
predictions of the currently favored cosmological model, namely Lambda-Cold-Dark-Matter (ACDM). 
We start by presenting the theoretical framework for such a comparison. It was flrst derived in 
the context of then-popular models like cold-d a rk-ma tter and isocurvature ba ryon, as described in 



detail in two classic pa pers: 



approach was taken by 



Juszkiewicz et al. 



Lavaux et al 



fll990h and 



Lahav et al. 



(119901 ). More recently, this 



(120101 ). who reconstructed the local peculiar velocity fleld (up 
to ~ 150 Mpc/Zi), applying the data from the 2MRS. The basic quantity for these comparisons is the 
joint probability distribution function for v and g, assumed to be a multivariate Gaussian. 

In our case, we want to flnd the expectation value for the amplitude of the acceleration of 
the Local Group knowing its peculiar velocity. The latter equals to fcMB = 622 ± 35km/s in the 
direction (/,6) = (272° ± 3°, 28° ± 5°). In ou r analysis we do not us e the directional constraint: 



the misalignment angle is integrated out (cf. 



Juszkiewicz et al 



1990l ). In the forthcoming paper. 



concerning the maximum-likelihood determination of P from 2MASS XSC, we will use the fact that 
the misalignment angle is known. Here, the theoretical prediction gives us the conditional velocity, 
Vc (in units of km/s), which is related to the acceleration via the l inear-theory relati on (ITT]) . The 



relevant formula for its amplitude has the form given by eq. (8a) of 
typo therein): 



Juszkiewicz et al. 



(1199 



o|) (note 



Vc = {Vm\vcMB) = 

1 - p2 + p^u^ 



erf 



pu 



pu [[2(l-p2)]l/2 

where we have used the following quantities: 



— (1 — p2) exp 

TT 



p^u^ 



2(1 -p2) 



Vm is the velocity induced by a single realization of the density fleld given the assumed power 
spectrum 



Vrr,. = 



f S{r)^W{r)d'r, (19) 
IT J 



4- 

with 5{r) and W{r) being respectively the density fleld and the window of the survey. The 
latter is included to mimic the velocity measured from a given survey; 
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0"^ is the predicted (ensemble average) variance of the velocity measured from the survey, given 
by 



oo 



^1= 'i T' Pik)w\k)dk, (20) 





with P{k) = (|5fcP) being the power spectrum of density fluctuations, k the wavevector and 
w{k) the observational window in Fourier space; 

U = Vm/o-m ; 

p is the correlation coefficient of Vm and fcMB 



oo 

2 



6 TT^ cr a 

where a is given by Eq. fl2U]) with ii;(A;) = 1 



/■ P(k)wik)dk, (21) 





The second moment of the conditional velocity is given by eq. (8b) of i Juszkiewicz et al.l (119901 ): 



/^c = (^^^I^^cmb) = 3a^(l - p^) + p^mV^ , (22) 

hence the variance is equal to 

ctI=Pc-v^,. (23) 

The presence of the factor /(fim) in Eqs. f|T9|) - fl20|) and the fact that we also have to include biasing b 
mean that in order to properly compare the conditional velocity Vc given by Eq. f|T8l) with the dipole 
d measured from the data, Eq. f[T^ . the latter has to be rescaled, using a best-fit /3 parameter, as in 
Eq. (fTTj) . In fact, this rescaling enables us to estimate /3 in a straightforward way; we come back to 
this point a few paragraphs later. 



6.2. Observational window 

In this subsection we discuss the proper observational window for our measurement. We start 
by noting that when calculating the peculiar gravitational acceleration of the Local Group from the 
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dipole of the all-sky galaxy distribution, in general there are two s chemes to postulate a relation 



Erdogdu et al. 



20061). The first 



between this distribution and that of the underlying mass (see e.g. 
scheme, called number weighting, assumes that the mass is distributed in the Universe as a continuous 
density field, which is sampled by galaxies in a Poissonian way. The second prescription, namely flux 
weighting, uses the assumption that all the mass in the Universe is locked to the mass of the halos of 
luminous galaxies. The 2MASS dipole as calculated in Sec. IHis of the fiux-weighted type. 

The observational window W{r) of a fiux-limited survey is just its selection function ip{r). In 
the number-weighted scheme, this function measures what fraction of galaxies located in the distance 
interval (r, r + dr) are included in the survey: 



oo 

/$(L)dL 



(24) 



where $(L) is the lumi nosity function and ^min is the limiting flux of the survey (e.g. 



1982 



Yahil et al 



Davis fc Huchra 



199ll ). In our case of the flux- weighted dipole, as given by Eq. f|T2|) . the appropriate 



window is the flux- weighted selection function, \l/(r), deflned in such a wa y that 1 — \l/(r) is t 



percentage of light frorn a distance r which is not visible for the survey (lErdogdu et al. 



he 



2006 



Chodorowski et al. 



20081 ). It is given by 



■qjfr) 



r L^L)dL 

^miii ^ ' 



(25) 



/-L$(L)dL ' 

where L^^in = 47rr^S'min- The window Ws{r) = \I^(r), where the subscript S emphasizes the 
dependence of the window on the flux limit of the catalo g, smoothly decreases to z ero with increased 
distance from the observer. For a detailed discussion see 



Chodorowski et al. 



(l2008f ). 



On the other hand, 'ideal' surveys are distance- (or volume-) limited: in this case there is no 
loss of signal up to the limiting distance R. Then the window function is simply unity for r < R 
and zero otherwise. It ha s thus a form of a Heav iside step function; for a spherical survey we have 



WR{r) = Qu{R-r) (e.g. 



Juszkiewicz et al 



19901 ). In order to better reconstruct the dipole of the 
galaxy distribution from a flux-limited catalog, volume-weighting of the survey is commonly mimicked 
by weighting individual gal axies by the inverse of the selection function, at the expense of inc reasing 



shot noise from large scales (IStrauss et al. 



19921 ) and at a price of a possibly large rocket effect (IKaiser 
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19871 ). The top hat is then the relevant window. In our case however, as we are deahng with angular 
(photometric) data only, we cannot estimate the selection function even in redshift space. Therefore, 
we do not weight galaxies in the sample (or in other words, we assign unit weights to them) and the 
relevant window function is given by (1251) . 

Th e Fourier-space counterpart of any observational window, wjk), is obtained from the real-space 



one as (iKaiser Sz Lahav 



1989 



Lahav et al. 



1990 



Juszkiewicz et al. 



19901 ) 



POO 

w{k) = k I W {r) iiikr) dr . 
Jo 



(26) 



where ji is the spherical Bessel function of the first kind of order 1. Note that this window is not the 
Fouri er transform of W{r). In our case of W{r) given by fl25l) . the Fourier-space window reads (eq. 



59 of 



Chodorowski et al. 



20081 ) 



Ws{k) = 1 



87rS'mir 



oo 

j sin(A;r)Linin'^'(^min)dr , 



(27) 



fc/„°°L$(L)dL 

u 

where S^i^ and Lmin are as above. The window Wg is a function of both the wavenumber k and the 
minimum flux Sram (hence the subscript s). On the other hand, it does not depend on the distance, 
as the latter is integrated out, together with absolute luminosities (with the use of the luminosity 
function). This means that when comparing observations to theoretical expectations, we never use 
the effective distance (1171) for the calculations of the conditional velocity Vc ( |T8l) and its variance 
cTc ( l23l) . Our observable is the minimum flux of the sample, related to maximum magnitude via 
'S'min = 'S'o lO^^-^^™^". For clarity however, in the relevant plots we prefer to show the results in terms 
of a linear scaling in the effective distance. 

In Figure [5] we plot the 2MASS flux-limited window with iTmax = 13.5 mag and compare it 
with two top-hat windows for distance-limited surveys: one with -Rmax = 296 Mpc/h and the other 
with -Rmax = 105 Mpc/h, which are respectively the effective distances of galaxies at the limit of 
the 2MASS and 2MRS Kma^ = 11.25 mag sample. The Fourier form of the top-hat window is 



( jJuszkiewicz et al. 



1990|) 



^2;th(A;) = l-jo(A;i?max) (28) 

with jo = sinx/x being the spherical Bessel function of the first kind of order 0. The oscillating 
behavior of this window in the Fourier space is a result of the rapid truncation in the real space. As we 
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k [/i/Mpc] 



Fig. 5. — Comparison of three observational windows in the Fourier space: the one of the 2MASS flux- 
hmited survey (thick red hne) for a maximum magnitude of 13.5 mag and two for distance-hmited 
surveys (top-hat) with -Rmax = 296 Mpc//i (blue thin line) and -Rmax = 105 Mpc//i (black dashed 
line). The former radius is the effective distance of galaxies at the limit of the 2MASS XSC; the latter 
corresponds to K = 11.25 mag, the limit of the 2MRS. 
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can see, the top-hat window with -Rmax = 296 Mpc/h, which asymptotically behaves in the same way 
as the 2MASS one with -ft'max = 13.5 mag, passes clearly more large-scale signal {k < 0.02 h/Mpc). 
This will result in much faster convergence of the dipole measured through the top-hat window than 
through the 2MASS flux- weighted one (see below). On the other hand, the top-hat window with 
a cut-off equivalent to the limit of 2MRS at 105 Mpc/h blocks almost all the signal already for 
k < 0.01 h/Mpc, as expected. 



6.3. Comparison with theoretical predictions 



We start by reminding that our whole analysis is performed within the linear theory. For that 
reason, we do not include non-linear effects in the correlation coefficient p (Eq. fUt) . Mathematically, 
this means that we set to unity two functions: the ratio of power spectra of density and velocity 



divergence, TZ(k), and th e cohe rence function, C{k). For n on-linear ana 



Cieciel^g &: Chodorowskil (120041 ). The latter, introduced by 



Strauss et al. 



y sis o: 



the former, see 



fll992l ) (originally called 



the decoherence function), can be interpreted as the cros s-correlation coefficient of Fourie r modes of 



gravity and velocity fields and was discussed in detail by 



Chodorowski fc Ciecie^d (120021 ). 



Our goal here is to compare the observational data to the expectations of the ACDM concordance 
model. Hence the power spectrum u sed in Egs. (I20l) and (|2T|) is the spectrum of cold dark matter, with 



baryon effects included, as given by 



Eisenstein&Hu 



(119981 ). In this framework, P{k) oc k'^'T^^k), 



where rig is the spectral index of primordial fluctuations and the transfer function T{k) depends on 
the parameters h, Q^, fib (density parameter of baryons) and (present value of root-mean-square 
density contrast of matter fluctuations within a sphere of 8 Mpc/h). In our c alculations we use 



201ll ): 



the following set of parameters obtained from WMAP seven-year observations (jLarson et al 
h = 0.71, n^h^ = 0.1335, Q^h^ = 0.02258, = 0.963 and as = 0.801. 

Before applying the observational data, in Figure [6] we compare the expected growth of the 
conditional velocity for the 2MASS window and for the top-hat case, both calculated for the ACDM 
model parameters. As we can see, for such a fiux-limited survey as 2MASS, the predicted velocity 
is far from converging to the limit of fcMB = 622km/s (horizontal line) even for r^s — 350 Mpc/h. 
In the top-hat case, the expected convergence is much faster, as could have been deduced already 
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Fig. 6. — Theoretically predicted growth of the conditional velocity of the Local Group for known f cmb 
using two different observational windows: the one of 2MASS (blue dotted line) and top-hat (green 
dashed line), both calculated for the ACDM as the underlying cosmological model. The horizontal 
solid line is the observed velocity of the Local Group with respect to the CMB. For comparison, the 
prediction for the currently disfavored Standard CDM model with the top-hat window is shown with 
the black dot-dashed line. 



- 31 - 



from the comparison of the windows, presented in Fig. [51 However, it should be noted that even 
for all-sky catalogs that include redshifts, for which we can effectively model the dataset with a 
top-hat window (like in the case of th e 2MRS), the c o nverg ence of the dipole is not likely before some 



200 Mpc//i, opposite to the results of lErdogdu et al.l (120061 ) . where it is claimed that the contribution 



from structure beyond 6000 km/s(= 60 Mpc/h) is negligible. Figure shows that for the latter 
distance, the conditional velocity for the top-hat window has reached less than 75% of its final value. 
On the other hand, for the 2MASS window, Vc attains 0.95 fcMB no sooner than for r ~ 470 Mpc/h 
(~ 14.5 mag in the Ks band), far beyond the completeness of the 2MASS XSC As an extreme case, 
in Figure Owe additionally plot the prediction for the once-popular 'Standard CDM' model, currently 
strongly disfavored by observations {h = 0.5, Q^n = 1, = 0.05, = 1, cxg = 1), with the top-hat 
window. Note how fast convergence would be expected in this case. 

The most important result of this Paper is presented in Figure [71 It shows the observed growth 
of the 2MASS dipole compared with the conditional velocity calculated from Eq. f[T8|) . As was already 
discussed at the end of Subsec. 16.11 we rescale observational data by a best-fit parameter f3 = VL^^ /h. 
The errorbars of the measured dipole were obtained from 1000-fold bootstrap resampling of the 
catalog. The /3 parameter was fitted by a minimum-x^ procedure, including both observational errors 
(from bootstrapping) and the theoretical (cosmic) variance. Note that formally we are 'over-fitting', 
as the datapoints are correlated with each other, due to an integral nature of the measured dipole. 

The observed growth, once rescaled, is well within la range of the theoretical prediction 
(apart from the datapoints for the smallest distances, where the number of galaxies is very small, 
hence the measurement is noisy). Remarkably, the best agreement is obtained in the range 
100 Mpc//i < Teff < 300 Mpc//i, i.e. between the extent of the 2MRS catalog and the 2MASS XSC 
completeness limit. The result of the fitting gives /3 = 0.38 with la confidence intervals of ±0.02. 
Owing to the considerations of Sec. [3l concerning the scatter introduced by the offset added when 
pa ssing from isophotal fiducial e lliptical aperture magnitudes to total fluxes, as well as due to findings 



of 



Bilicki &: Chodorowskil (|2010| ) on the bias introduced by the fact of possibly improper masking of 



the Local Void, we double this formal error. Hence our estimate is 



(3 = 0.38 ±0.04 



(29) 
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Fig. 7. — Growth of the 2MASS clustering dipole (red points with la errorbars), compared with the 
theoretical expectation for conditional LG velocity in the ACDM model (solid blue line with dotted 
blue lines showing Icr variances). The datapoints were reseated by the /3 parameter (see text). The 
horizontal solid black line is the Local Group velocity with respect to the CMB. Vertical lines give the 
limits of the 2MRS and 2MASS XSC catalogs. Note the general consistency of the observations with 
theoretical expectations. 
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This value is in agreement with the analysis of 



Erdogdu et al. 



(120061 ) ■ where also data from 2MASS 



were used, although from much smaller scales. The clustering dipole was calculated there from 
a subsample of the 2MASS XSC with measured redshifts, namely 2MASS Redshift Survey, then 
complete up to ~ 100 Mpc//i (maximum magnitude in the Kg band equal to 11.25 mag). Their result. 



based on claimed conv ergence of the fl u x-wei ghted dipole, was /3 = 0.40 ± 0.09. A somewhat larger 

2^ 



value was obtained by iPike fc HudsonI (120051 ). where the peculiar velocity field within 65 Mpc//i 
predicted from 2MASS photometry and public redshift data was compared to three independent 
peculiar velocity surveys based on Type la supernovae, surface brightness fluctuations in elliptical 
galaxies, and Tully-Fisher distances to spiral galaxies. The best-fit from this comparison was 
(3 = 0.49 ± 0-0 4 - Oil the other hand, our value of /3 agrees within the errors with that obtained 



by 



Davis et al. 



(120111 ). who reconstructed the cosmological large scale flows in the nearby Universe 
using the SFI++ sample of Tully-Fisher measurements of galaxies and compared it with the whole 
sky distribution o f galaxies in the 2MR S to derive 0.28 < (3 < 0.37 (68.3% confidence). Finally, 



in a recent paper. 



Nusser et al. 



(l2011al ) used the galaxy distribution in the 2MRS to solve for the 



peculiar velocity field, estimated absolute magnitudes of galaxies and constrained /3 by minimizing 
their spread. They found a value in excellent agreement with our determination: /3 = 0.35 ± 0.1. 



7. Summary and conclusions 

The motion of the Local Group of galaxies through the Universe can be used as a tool to 
constrain cosmological parameters. One of the applications is to estimate the peculiar acceleration 
of the LG and compare it within the linear theory with its peculiar velocity, known from the dipole 
anisotropy of the cosmic microwave background. This comparison may be used to directly measure 
/3 = Q'^^^/b, and hence the cosmological parameter of matter density, Q^- 

In this paper, as the estimator of the LG acceleration we have used the clustering dipole of 
the galaxies from the 2MASS Extended Source Catalog, which contains positions and fluxes of 
almost one million sources. As we wanted to reach as far from the Local Group as possible with 
the state-of-the-art observations, we have decided to trade the advantages of redshift surveys for the 
huge number of objects and unprecedented sky coverage of the 2MASS XSC. This allowed us to 
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measure the dipole up to 300 Mpc/Zi, which is almost 3 times farther than for the deepest existing 
all-sky survey of galaxies with measured redshifts, the 2MASS Redshift Survey (2MRS, a subset of 
the 2MASS XSC, containing ~ 43,000 galaxies). The price to pay was our inability to weight galaxies 
with the inverse of their selection function (which would require redshifts to be known), so we could 
not directly mimic an 'ideal', distance-limited survey, with our data. On the other hand, the flux 
dipole does not suffer from redshift-space distortions, and in particular from Kaiser's rocket effect. 

We have focused mainly on the issue of convergence of the dipole, which has been a subject of 
debate for more than two decades now. We have shown that the flux dipole of the 2MASS XSC 
does not converge up to the completeness limit of the sample, which is 13.5 mag in the Ks band. 
Moreover, we find that beyond an effective distance of some 150 Mpc//i, the growth is induced mainly 
in the negative Galactic Cartesian y direction. We have checked that this behavior is not due to 
systematic effects related to masking out and artificially filling the Zone of Avoidance. Additionally, 
the misalignment angle between the measured dipole and the velocity of the LG is found to be of the 
order of 20°. We also find that the misalignment comes mainly from some nearby galaxies, bright 

could optimize the analysis, as was already 



Chodorowski et al. 



(120081 ). We plan to come back to this issue in the 



in the near infrared. This is a hint that their remova 
theoretically investigated by 
future. 

Our inability to weight galaxies and to mimic a distance-limited survey and the lack of 
convergence of the dipole do not allow us to measure the parameter /3 from direct vlg ^ Qlg 



comparison. However, by comparing the observed growth of the clustering dipo^ 



expectations for known fcMB? in the framework of 



Juszkiewicz et al. 



fll990h and 



e with theoretical 



Lahav et al 



fll990h 



we find consistency with the predictions of the ACDM model, once the proper observational window 
of the 2MASS fiux-limited catalog has been included. This fact enables us to measure the /3 parameter 
from this comparison alone and the result is 



(5 = 0.38 ±0.04 . 



(30) 



This value is in agreement with and has a better prec ision than an e a rlier 



used data from 2MASS, although from smaller scales: 



Erdogdu et al. 



analysis, which equally 



fl2006h found (3 = 0.40 ± 0.09 



from the clustering dipole of the 2MASS Redshift Survey, with a depth of ~ 100 Mpc/h. 



- 35 - 



Having estimated the value of /3, we can go further and d erive the dens i ty pa rameter fim provided 



Mailer et al 



(l2005l ). who calculated the 



that we know the linear bias, b. The latter was estimated by 
angular correlation function of galaxies in 2MASS and inverted it using singular value decomposition 
to measure the three-dimensional power spectrum. A fit of CDM-type power spectra in the linear 
regime allowed them to constrain the i^<j-band linear bias as = 1.1 ± 0.2. Using this value, we 
obtain the following estimate of the density parameter: 



0.20 ±0.08 



(31) 



This result, supporting the fact that the density of non-relativistic matter in the Universe is well 



bellow t 
that by 



le critical one, is con sistent with other independent dynamical estimates of Q^, such as e.g. 



Feldman et al. 



(120031 ). who used mean relative peculiar velocity measurements for pairs of 



galaxies and obtained Qrn = 0.30 



h0.17 
-0.07- 



Our findings can be verified in at least two ways. The first would be to examine the growth of the 
2MASS dipole with redshifts of galaxies as their distance estimates, using proper weights. In the near 
future we cannot however hope for a uniform and deep enough sample of spectroscopic redshifts for 
the catal og, even if the 2MA SS Redshift Survey is continued (as it was planned to reach Kg = 12.25 



mag, see 



Huchra et al. 



20051 ). The calculation may be thus feasible only for photo-^'s, if they are 



available for the whole 2MASS XSC A promising direction towards this goal may be to cross-cor relate 



2MASS with the data from the Wide-field Infrare d Survey Explorer (WISE, iWright et al. 



are currently partially available (jCutri et al. 



2010h that 



20111 ) and are expected to be fully released in 2012. The 



other possible verification method could be for instance to examine the behavior of the differential 
dipole, i.e. to compare the increments of the growth with theoretical predictions of the ACDM model. 
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A. Effective distance for a given flux 

In this Appendix we show how we calculate the effective distance of galaxies with a given flux 5", 
knowing their luminosity function (LF) in the given band, ^{L), but not knowing their redshifts. Let 
us start by deriving the effective depth of a fiux-limited sample. 

The number of galaxies cLA^ in a volume element dV = 4 vr dr of a spherical sample limited by 
minimum flux 6*111111 (equivalent in the K band to some limiting magnitude -ft'max) is given by 

r+00 

dA/" = 4 TT rMr / $(L)dL, (Al) 

where L^i^ = 4 minimum luminosity of galaxies in the sample at a given distance r. 

The mean depth of the sample is defined as 



which, for the Schechter LF (jSchechter 
L*, gives 



f+00 ^ <w ^ 

Jo dr ^' 



19761 ) with a faint-end slope a and characteristic luminosity 



(R) = ^00 .::: .r":?:: > (A3) 



/+°°r3r(l + a,4 7rr^ W-^*)dr 
/+°"r2r(l + a,4 7rr2 5iiiin/L,)dr 
where r{a,x) is the upper incomplete Gamma function. 



Using the distribution (lAip . we can also easily derive the median depth of a fiux-limited sample. 
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which we shall denote as R. We obtain it by solving the integral equation: 



R 



dr 



+00 



dr 



dr 



dr . 



(A4) 



R 



The following Table [T] presents the mean and median dep 
and 2MASS XSC (K^ax = 13.5 mag) for the LF as given by 



th of the 2MRS (K , 



Jones et al. 



11.25 mag) 



(120061 ). The difference 



between the mean and the median is in that case of the order of a few percent. 



survey 


mean depth (R) 


median depth R 


2MRS 


82 Mpc/h 


76 Mpc/h 


2MASS XSC 


232 Mpc/h 


215 Mpc/h 



Table 1: Mean and median effective depth of flux-limited surveys: the 2MASS Redshift Survey (2MRS, 
Kma.x = 11.25 mag) and the 2MASS Extended Source Catalog (XSC, i^max = 13.5 mag), calculated 



from Eqs. ( 1A3I) and ( ]A4[) . respectively, 
a = -1.16 and = -23.83 + 5 log h, 



or the /T-band lum inosity function in the Schechter form with 



Jones et al. 



(120061). 



The effective depth of a flux-limited survey is however not a good measure of the effective 
distance of a thin shell in 'flux space' (in which the third coordinate is the flux, by analogy with the 
redshift space) that we seek. What we need is the mean value of distances of all galaxies with given 
flux S. The mean that we calculate will be thus a conditional one. This derivation is qualit atively 



the sa me as for the mean redshift of galaxies with a given flux, presented on pages 120-121 of 



Peebles 



(|l993h . 



We start by deriving the joint probability distribution of galaxy distances r and fluxes S. It 
is easily obtained by differentiating Eq. ( lAll) with respect to luminosity L. Hence, the differential 
number of galaxies with a LF $(-L) in a volume element 6V is given by 



S^N = ^(L) 5LSV 

with 6V = ATcr"^ 6r. Now, passing from luminosity to flux, L = Anr^S, for fixed r we have 
SL = 47rr^ 6S. This gives the joint probability as 



(A5) 



p(r, S) 



dr dS 



16 TT^ r^$(47rr^5). 



(A6) 
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The conditional probability for r given 5* is 



VMS) = ^ = ^6vr^r^$(W^) 

' ' piS) /+~16 7r2r4<|.(47rr25)dr ' ^ ' 

hence the conditional mean for r given S will be 

Jq $(47rr2S^) dr 

For the Schechter form of the luminosity function, let us define a characteristic distance rs, which is 
a distance to a galaxy with a given flux S and a characteristic luminosity L^: 



Then for the Schechter LF, Eq. (lASP simplifies to 



r 5 = rc ^ ; AlO 
r(Q; + 5/2) 



with r(a) being the Gamma function. 



Similarly, the conditional median for r given 5, denoted here as rs, will be obtained from the 
implicit equation 



rs +0O 



r*<l>(47rr^5)dr = y $(47rr^5) dr , (All) 

fs 

which for the Schechter LF is equivalent to solving 

with respect to fs, where 7(0,0;) and T{a, x) are respectively the lower and upper incomplete Gamma 
functions. 



In the particular case of the K band, taking a = —1.16 ( 1 Jones et al.ll2006l ). we obtain 

{r)s = 1.056 rs (A13) 

and 

fs = lM3rs. (A14) 

Note that especially the median conditional distance gives a value very close to the characteristic 
distance rs, which could be the 'first-guess' effective distance (Tully 2008, private communication). 
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For the purpose of the present paper, we have decided to use the conditional median as our 
measure of the effective distance of galaxies with a given flux. This particular choice is somewhat 
arbitrary and influences scalings of top axes in Figs. [2],|H[6]&[71 However, since the difference between 
the conditional mean and median is very small (about 5%), this choice has a negligible impact on 
general conclusions and results of this paper. 

Relating the flux to the K magnitude via Eq. f|T4l) . S = Sq lO^"'"^^, and owing to the 
luminosity-magnitude relation L* = 47r(10 pc)^5'o lO^^ '^*-'^*, we obtain 



reff = fs = 1.013 X 10°-2(^-^*)-^Mpc 



(A15) 



so finally, for M, = -23.83 + 51og/i fl Jones et al.ll2006h 



Tcff = 0.591 X 10°-^^ Mpc//i 



(A16) 



For example, galaxies at the limit of the 2MASS catalog (with K = 13.5 mag) are assigned an 
effective distance of Ves = 296 Mpc/h. 
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